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We report low temperature magnetotransport measurements on a high mobility (/x = 
325 000 cm 2 /V sec) 2D electron system on a H-terminated Si(lll) surface. While low magnetic 
field data indicate a six-fold valley degenerate system, we observe the integral quantum Hall effect 
at all filling factors v < 6 and find that v = 2 develops in an unusually narrow temperature range. 
An extended, exclusively even numerator, fractional quantum Hall hierarchy occurs surrounding 
v — 3/2, consistent with two-fold valley-degenerate composite fermions (CFs). We determine acti- 
vation energies and estimate the CF mass. 

PACS numbers: 73.40.-c, 73.43.-f, 71.70.Gm 



Multicomponent systems in a single Landau level have 
generated interest due to the possibilities for novel corre- 
lated ground states in the integer and fractional quantum 
Hall (FQH) regimes when the energies of the component 
states become degenerate. Early experiments focused on 
measurements of engineered GaAs materials where the 
spin splitting could be reduced to zero [1, 2j. Subsequent 
development of AlAs quantum wells with a tunable valley 
degeneracy allowed the study of a spin-like degeneracy in 
the same limit [•>] . Recently, there has been great interest 
in the sublattice (valley) degeneracy in graphene; exper- 
iments show that the valley symmetry affects the FQH 
hierarchy [4, 5] and that valley ferromagnetism occurs 
when one of two degenerate valleys is occupied [G] . 

Measurements on silicon, the first multi- valley system 
to be considered theoretically [8, 9], had been hampered 
by high disorder. Lately, however, Si(100)/SiGe het- 
erostructures have shown mobilities up to ~ 10 6 cm 2 /Vs 
[10]. Nonetheless, Si(100) is known to have an intrin- 
sic valley splitting due to surface effects [11—14]. The 
case of 2D electrons on (111) oriented silicon surfaces, 
which have three pairs of opposite momentum (fc) val- 
leys, is especially interesting. As opposed to either AlAs 
or Si (100), valley pairs with ±fc symmetry in Si(lll) can- 
not be broken within the effective mass approximation or 
by a surface potential [!)], similar to the case of valleys 
in graphene. Additionally, both AlAs and Si(lll) ex- 
hibit anisotropic mass tensors; in AlAs, this anisotropy 
arguably transfers to composite fermions [ 1 5] . Novel bro- 
ken symmetry states have been predicted at integer filling 
factor, v, when the Fermi energy, Ep. lies between two 
valleys with different mass tensors [16]. 

In this Letter we present transport data on a very high 
mobility electron system (fi = 325 000 cm 2 /Vs at a tem- 
perature T = 90 mK and density n s = 4.15 x 10 11 cm~ 2 ) 
on a hydrogen-terminated Si(lll) surface. We discuss 
the anisotropy between orthogonal current directions and 
observe an extended FQH hierarchy around v = 3/2. 
We argue that the FQH hierarchy is consistent with the 
SU(2) symmetry of a valley-degenerate ground state, and 



estimate the mass of composite fermions. Finally, we 
present evidence for many-body interactions affecting ac- 
tivation energies of several integer filling factors. We find 
that the development of v = 2 occurs in an unusually nar- 
row temperature range, and thus may signal a transition 
to broken symmetry valley states. 

For our samples, we replace the Si-SiC>2 interface of 
a typical MOSFET with an interface between H termi- 
nated Si(lll) and vacuum in order to remove the effects 
of strain and dangling bonds from the Si(lll) surface 
[17]. This is done with the help of an encapsulating 
silicon-on-insulator piece which is bonded via van der 
Waals forces to a high resistivity H-Si(lll) substrate and 
forms a gate. The sample studied in the present work is 
fabricated in the same way as previous ones, with fur- 
ther optimized cleaning and annealing steps [18]. Four- 
terminal resistance in a square van der Pauw geometry 
is denned as Rij,i m = Vi m /hji with R xx = #12,34 and 
R yy = #24,13 oriented along the [110] and [112] direc- 
tions respectively (see inset to Fig. lb). Hall traces, p xy , 
are the averages of #14,23 and #23.14 to prevent mixing. 
The data were obtained at densities between n s = 3.7 
and 5.7 x 10 n cm~ 2 , adjustable via a gate voltage, in 
both a 3 He system and 3 He/ 4 He dilution system with 
base temperatures of 280 mK and 90 mK, respectively. 
The density range is limited due to nonlinear contact re- 
sistances at low n s and gate leakage at high n s . Magneto- 
transport measurements were performed using standard 
lock-in techniques with a typical excitation of 1 — 25 nA 
at 5 Hz in a magnetic field up to 12 T. 

Figure la shows R xx , R yy and p xy versus magnetic field 
B taken at T = 90 mK and density n s = 3.75 x 10 11 cm' 2 . 
Shubnikov-de Haas oscillations (SdHO) are visible down 
to about B ~ 0.15 T. At low B and before the onset 
of spin splitting, minima occur at Ais xx — 11.65 ± 0.04, 
while Ai/yy — 12.6 ±0.1, roughly consistent with a valley 
degeneracy g v = 6 and a spin degeneracy g s = 2. Due 
to the anisotropy between R xx and R yy (discussed be- 
low), we do not expect perfect agreement with Av = 12. 
By B sa 1.3 T [y — 12), only the lowest Landau level 



2 




B(T) 

Figure 1. (a) R xx , R yy (left), and p xy (right) vs. B at 
T = 90 mK and n s = 3.75 x 10 11 cm -2 with filling factors 
listed on the top axis; the fractional states at u = |, and 
| have plateaus in p xy . (b) Calculated anisotropy p xx /p yy for 
the same data as (a). Dotted lines indicate regions near fill- 
ing factor minima, where R is small and leads to large errors. 
The dot-dashed lines show predicted anisotropy if only a sin- 
gle valley pair contributes to conduction (see inset and text). 
Inset: Schematic of the sample contacts and valleys oriented 
along major orthogonal directions of our Si(lll) samples. 



is occupied. As B increases, the six-fold valley degener- 
acy breaks. Below v — 8 all integer filling factors have 
minima, with the even states (y — 6, 4, and 2) being 
stronger than the odd (y — 3, 5 and 7). Below v = 2, we 
observe many fractional states with strong minima and 
associated plateaus at v = 8/5, 10/7, and 4/3. 

While R xx sa R yy at B — 0, strong anisotropy appears 
for B > 0. Such anisotropy is to be expected in Si(lll) 
when valleys with different mass tensors have an unequal 
density of states (DOS) at Ep. If only a single valley is 
being occupied, and scattering is isotropic, then the the- 
ory of Ando [19] predicts that for m x > m y , p X x/Pyy = 
m x /m y , where m x and m y are the principal masses of 
the valley, and the symmetry axes of the resistivity ten- 
sor are those of the valley. Using notation appropriate 
for our system, the model predicts a resistivity anisotropy 
that lies within the range 0.28 < p xx /p yy < 3.55 using 
m x = 0.190m e and m y = 0.674m e . 



As noted by Simon in the context of anisotropics at 
v — 9/2 [20], the resistance anisotropy can greatly ex- 
ceed the resistivity anisotropy in van der Pauw geometry 
measurements if the symmetry axes of the resistivity ten- 
sor are aligned with the device edges [21]. This effect will 
be largest in our samples when the "C" valleys are being 
occupied (see the inset to Fig. lb). The largest resis- 
tance anisotropies we observe in our data occur in the 
range 4 < v < 6 where R yy /R xx w 95. To account for 
this we assume that the symmetry axes of p are those of 
the "C" valleys and calculate p from R using formulas of 
Versnel for a circular van der Pauw geometry [ ]. The 
calculated minimum ratio p xx / 1 p yy = 0.23 (see Figure lb) 
agrees well with the Ando result when a single "C" valley 
is being occupied. At lower B, the system is well de- 
scribed by isotropic conduction for 6 < v < 8, consistent 
with the approximation that six valleys are being equally 
occupied in this range. As B is reduced further, however, 
the anisotropy continues to oscillate due to SdHO. The 
relative valley occupation is unclear for 2 < v < 4 and 
therefore we cannot predict the anisotropy. 

When v < 2 we expect that only the final pair of val- 
leys are being occupied and observe a nearly constant 
resistance anisotropy. The ratio is smaller than that ob- 
served at lower B, and the value of R xx /R yy crosses from 
< 1 to > 1 for v < 2.25 (B > 6.7 T). These observations 
can be explained if either only "A" or "B" valleys are be- 
ing occupied. Using the equations of Versnel [22], p is 
determined from R assuming that the symmetry axes of 
p are oriented 120° (or equivalently 60°) away from the 
sample edges. The result, plotted in Fig. lb again gives a 
maximum resistivity anisotropy that agrees well with the 
Ando prediction for a single valley [19]. The anisotropy 
present in our square samples can therefore be related to 
the underlying anisotropic mass tensor combined with an 
unequal valley DOS at Ep. 

We focus now on the high B regime of the fractional 
quantum Hall effect. Figure 2 shows the temperature 
dependence of R xx in the FQH regime below v = 2. In 
this range, we observe 9 states: v = 8/5, 14/9, 20/13, 
22/15, 16/11, 10/7, 4/3, as well as 6/5 and 14/11 (not 
shown). The hierarchy of observed states is exactly that 
predicted by an SU(2) symmetry in which the two-fold 
valley degeneracy of electrons is preserved, leading to the 
creation of composite fermions (CF) with two and four 
attached vortices ( 2 CF and 4 CF, respectively) [23]. From 
the anisotropy analysis above, the two-fold degeneracy 
stems from the occupancy of either an "A" or "B" valley 
pair. For states below v = 1, the filling factor v* of 
CFs is given by v — v* /(2pv* ±1) (where p = 1 and 
2 for 2 CF and 4 CF, respectively). On the other hand, 
the hierarchy of fractional states around v = 3/2, which 
are the hole-symmetric equivalents to the v = 1/2 states, 
is given by v = 2 - v*/(2pv* ± 1) [23]. With an SU(2) 
symmetry, only v* — 2, 4, 6, . . . are expected, giving rise 
to the hierarchy visible in Fig. 2a. A similar, though 
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Figure 2. (a) R xx vs. B for n s = 3.75x 10 11 cm -2 . Seven FQH 
states are visible at T = 110 mK; the filling factors are labeled 
on the top axis. Top panel: p xy versus B for T — 110 mK. 
The derivative is also shown to accentuate the fine structure, 
matching the R xx minima, (b) Mass of composite fermions 
assuming that the gaps probed by Shubnikov-de Haas (SdH) 
oscillations have the form Acf = huicF- Error bars are the 
95% confidence bars from fits to SdH data. On the right 
axis, the gap energy is plotted for each mass. Open symbols 
indicate values extracted from maxima. The two lines are fits 
to available gaps; the intercept for the fits is about — 0.09K. 
The gap value at v — 4/3 is calculated from activation energy 
data, and is consistent with the extrapolation of the mass. 



smaller, hierarchy was observed recently in graphene with 
the same conclusion [ ]. Finally, v — 6/5 and 14/11 
relate to 4 CFs via v* =4/3 and 8/5, respectively, using 
2 — v* /{2v* — 1); the 4 CF hierarchy simply reflects the 
FQH states of 2 CFs. We note that we observe no evidence 
for a v = 5/3 state, similar to recent experiments in both 
graphene [5] and Si/SiGe [10]. However, the 5/3 state 
is visible in nominally doubly-degenerate AlAs, probably 
due to local strains [16, 24, 25]. 

To estimate the mass of the composite fermions, we use 
the theory of SdHO and apply it to v = 3/2 by trans- 
forming to an effective magnetic field B e g = 3 (B — B 3 / 2 ) 
where B 3 / 2 — 10.34 T (calculated from the density) [26]. 
Although we could use either activation energy measure- 
ments or SdHO to find the gaps, analysis based on SdHO 
takes into account intrinsic level broadening. We first 
compute the amplitude AR of the oscillations using a lin- 



ear interpolation of the minima and maxima (see [27]). 
We then use AR oc Rq exp {-27r 2 7VA}£/sinh£, where 
£ = 2tt 2 T/A, Rq is the (temperature dependent) resis- 
tance at B c ff — to find the gap A. The Dingle temper- 
ature To characterizes the intrinsic level broadening. 



We expressly assume that Acf 



CF 



,CF 



e\B e ff\ /mcF and extract mcF directly; the results are 
shown in Fig. 2b; m CF /m e = 0.98 ± 0.03 (B cS < 0), and 
mcF/me = 1-20 ± 0.03 (B c g > 0). We also plot the gap 
energy Acf as a function of B e g, and fit a line to the re- 
sulting energies. The low intercept of —0.09 ± 0.07 K for 
-B e ff is consistent with the calculated CF mass. The acti- 
vation energy at v = 4/3, measured from higher T data, 
is consistent with a constant value for mcF- Additionally, 
we find that all of the FQH data collapses onto a single 
line for Ai?sinh£/i? £ vs. l/B with Tg F = 180 ± 8mK. 
This is more than twice as large as the value for electrons, 
where Tjj = 78 ± 3 mK, determined at low B using the 
same approach. 

While our data is well modeled with a two-fold valley 
degeneracy, a total- independent valley splitting which 
polarizes composite fermions, similar to effects due to 
spin polarization [26], could explain the difference in 
CF mass between ±B c g. Such a valley splitting would 
be consistent with that found in other silicon devices 
[28, 29]. Additionally, this splitting, which is small and 
has an upper bound of To at B = T, would also ap- 
pear in the FQH hierarchy. The derivative in p xy in Fig. 2 
shows a slight minimum at approximately v = 7/5, which 
would be the strongest odd- numerator state for v < 3/2 
for 2 CFs. Due to the weakness of the state, however, 
we conclude that the valley splitting is less than Tq F at 
£> c n = 0. A significant splitting would also be observed at 
v > 2; a weak minimum is seen at v — 7/3, and is the only 
odd-numerator state visible. The strength of this state 
is nonetheless consistent with fractional states due to in- 
creasing particle (as opposed to hole) population; mea- 
surements on graphene similarly show a hierarchy of the 
1 /3 family of states [4] . Finally, we note that while our es- 
timates for the CF mass are not entirely consistent with 
recent measurements on both GaAs [26, 30] and ZnO 
heterostructures [31], there are several sample-specific ef- 
fects that could renormalize the CF mass, including the 
finite perpendicular extent, z±, of the electron gas and 
Landau level mixing [32, 33]. The modification of short- 
range interactions due to the large dielectric-constant, e, 
mismatch at the surface is particularly relevant for our 
devices. 

We turn now to gaps at integer filling factors; we use 
the T dependence of R xx and R yy to calculate the acti- 
vation energies via R oc exp (— A !/ /2fcsT). By changing 
n s , we are able to measure the gap A„ as a function of 
B: Figure 3a shows the results of such a calculation in 
units of the Coulomb energy Eq = e 2 /47ree ^B, where 
Ib = y/h/eB and we use e = (eg; + eyac) /2 = 6.25, 
which is appropriate for interaction ranges 3> z±_ . Due to 
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Figure 3. (a) Activation energy measurements as a function 
of B in units of Coulomb energy. Vertical bars show the 95 % 
confidence errors from fits to the linear portions of Arrhenius 
plots. Solid symbols indicate data from R yy measurements, 
while empty symbols are from R xx data. The v = 2 data 
shows an enhanced energy gap greater than the Landau level 
spacing and indicates the scatter between R xx and R yy data. 
The v — 6 data is used to extract an effective g* = 3.63 ± 
0.06. (b) Plot of the linear portion of a yy (representative of 
both directions) versus A„/T, showing the difference between 
the critical conductivities a yy at different filling factors and 
densities. For v = 2, a yy is between 10 and 16 orders of 
magnitude larger than the other filling factors (regardless of 
their density) , and the color of each data point is the same in 
parts (a) and (b). 



anisotropy and a lack of large resistance range in R XXl we 
show only A^ for most of the filling factors; for v = 2 we 
show both transport directions as an example of scatter 
in the gap energy. 

With the assumption that the valley splitting is smaller 
than the Zeeman gap, the v — 6 activation energy is 
interpreted as Ag = Ez = g*fJ-BBtot, from which we 
estimate g* = 3.63 ± 0.06, consistent with previous mea- 
surements [17, 18]. The presence of minima at all filling 
factors v < 6 suggests a independent valley splitting, 
similar to that observed in SiGe heterostructures [28, 29]. 
Indeed, the appearance of odd v < 6 shows that high 
magnetic fields increasingly break the two-fold degener- 



acy of opposite k valleys. The estimate at v = 5 is based 
on a "strength" (S) of the state defined as the ratio of 
the resistance minimum to the average of the adjacent 
maxima [10]. From the T dependence, we can estimate 
a quasi-gap; the state is weaker than any other filling 
factor v < 6. For v = 1, we can only estimate an upper 
and lower bound while noting that the minimum remains 
visible at T = IK (see Fig. 3a). The odd filling factor 
gaps, which are much greater than Tjy, indicate that the 
splitting is likely due to many-body effects. 

At v — 2 one pair of valleys is completely filled, so 
added electrons must occupy a new valley with a dif- 
ferent mass tensor. In this regime, Abanin et al. [16] 
predict that novel nematic phases are favored, with the 
sample broken up into domains of differing valley po- 
larization. Experimentally, the temperature dependent 
behavior at v = 2 is strikingly different from other fill- 
ing factors. Indeed, we observe it to be characterized 
by a very narrow T range where R is thermally acti- 
vated, leading to a large extrapolated R for 1/T — > in 
an Arrhenius plot of the data. In Figure 3b, we show 
the thermally activated portion of a yy as a function of 
A v /T for the integer filling factors. The critical con- 



a yy (1/T — >• 0) in the IQH regime is ide- 



ductivity Oy y 

ally ~ 2e 2 /h, with reductions possible due to short-range 
scattering and screening [34-36]. Compared to a yy at 
higher filling factors, the critical conductivity at v = 2 
is 10 to 16 orders of magnitude larger. Abanin et al. ar- 
gue that in an anisotropic valley-degenerate system, the 
transport mechanism in the limit of valley-polarized do- 
mains is variable-range hopping due to the low T localiza- 
tion of edge currents along domain walls separating areas 
of different valley polarization. Fits to different models, 
including a — (cr c */T) exp {A l/ /2fc B T}, reproduce the 
large energy gap and do not reduce the discrepancy in 
a yy P^]' To the best of our knowledge, the thermal be- 
havior reported here for v = 2 has not been seen in any 
other QHE state. 

In summary, we have shown evidence of electron- 
electron interactions in a high mobility Si(lll) system. 
The six-fold valley degeneracy breaks at high magnetic 
fields into an apparent SU(2) symmetry reflected by the 
fractional quantum Hall state hierarchy. We estimated 
the mass of composite fermions near v — 3/2 by assum- 
ing that the Shubnikov-de Haas oscillation gaps can be 
interpreted cyclotron energy. Further experiments 
are necessary. First, hexagonal samples would allow un- 
ambiguous measurements of p on Si(lll) surfaces. Sec- 
ond, higher magnetic fields are necessary to probe the 
v < 1 regime at similar densities, which would shed fur- 
ther light on the valley degeneracy of composite fermions. 
Finally, tilted magnetic fields would introduce controlled 
valley splitting of incquivalent valleys [ , ] and in par- 
ticular would help shed light on the behavior at v = 2. 
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